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Time-reversal symmetry is important to optics. In linear optics, a time-reversed process 
to laser emission enables total absorption of coherent light fields into an optical cavity of 
loss by time-reversing the original gain medium. In nonlinear optics, time symmetry exists 
for some well-known processes such as parametric up/down conversion, sum/difference fre- 
quency generation, however, combining exact time-reversal symmetry with nonlinear wave 
mixings is yet explored till now. Here, we demonstrate time reversed wave mixings for sec- 
ond harmonic generation (SHG) and optical parametric amplification (OPA). This enables 
us to observe the annihilation of coherent beams under time-reversal symmetry by varying 
the relative phase of the incident fields. Our study offers new avenues for flexible control 
in nonlinear optics and potential applications in efficient wavelength conversion, all-optical 
computing. 

Time-reversal symmetry dominates in many physical systems; it allows a physical process to 
reverse in a backward direction of time. This powerful symmetry enables many practical applica- 
tions under such reversible principle, with examples including spin reversing in nuclear magnetic 
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resonance (NMR) imaging^, acoustic and electromagnetic focusing using time-reversal mirroP^, 
phase conjugation mirror in opticd^E, etc. Especially in optics, light transmission in random 
medium can be enhanced by several orders by exploring the time-reversal symmetry simi- 
larly for enhancing second harmonic signals in random mediumP-R Also a quantum splitter can be 
obtained by time-reversed Hong-Ou-Mandel interferometer^. Recently, a novel concept named 
"coherent perfect absorber" (CPA) EfflS that explores time-reversed process to laser emission has 
shown that: incident coherent optical fields can be perfectly absorbed by a time-reversed optical 
cavity by replacing the gain with equal amount of loss. Also the incident fields and frequency 
should coincide with those of the corresponding lasing modes with gain under time symmetry. 
Since then, many efforts have been made to study CPA properties with different geometries 
However, most of them remain in linear optics regime. Nonlinear version of CPA has been the- 
oretically proposed to investigate signal and idler beams' phase varying dynamics in presentence 
of pumping beam under time reversed optical parametric oscillation (OPO) scheme These 
studies of time symmetry have been attracting increasing attention, since it provides alternative 
and substantial ways to manipulate light in nonlinear regime. Here, we experimentally exam the 
time-reversal symmetry for two of classical nonlinear wave mixing processes: SHG and OPA, 
characterize their nonlinear properties as opposite to their time-reversal counterparts, and reveal 
the nontrivial dynamics of phase varying in time-reversed nonlinear wave mixing schemes. 

In nonlinear optics, several pairs of wave mixing phenomena seem to be reversible in time. 
For example, sum-frequency generation (SFG) converts two photons into one, natively the opposite 
scenario, difference-frequency generation (DFG) that splits one to two seems to be the counterpart 
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to the former one in time. However, time symmetry requires the exact inversion of all physical 
parameters with respect to time. These pair processes have dramatic different initial conditions in 
term of pumping intensity, polarization, phase, etc. In order to achieve time-reversal, one has to 
consider them all in time. 

The SHG process converts a fundamental wave (FW) into its second harmonic (SH) through 
nonlinear response of a medium. However, the energy flow can oscillate between the fundamental 
and harmonic waves, determined by the phase mismatch condition and phase difference of the 
interacting waves. Under microscopic theme, this energy oscillation between FW and SH results 
from interference effects due to nonlinear polarizations of electrical dipoles. The rates of such up 
and down conversions are determined by their nonlinear susceptibilities, according to permutation 
symmetry \ 2l \ 22 \ they are equivalent. The rest conversion relies on the phase relationship. Below 
the depleting pump limit, FW and SH can be described by two coupled linear equations. Even 
reversed in time , this coupling still obeys the same physical system. For a non-seeded scheme, SH 
starts to grow from quantum fluctuations draining the energy from FW along the propagation. The 
exact time reversal of the problem can lead SH to none for the initial stage, similar to CPA in the 
linear case. 

In this letter, we exam time reversed processes for SHG and OPA. These backward nonlinear 
wave mixings allow the unique property of annihilation of coherent beams in a nonlinear quadratic 
medium by time reversal. Unlike the case of CPA in linear regime where incident fields are totally 
absorbed and converted into head 13 ! 16 !, here annihilation of incident fields can lead to the generation 
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of new fields, such backward parametric interactions may have a future for efficient wavelength 
conversion for better long-wavelength detection, e.g. mid-IR, THz. More interestingly, a flexible 
phase control can be achieved to probe the nonlinear dynamics during the wave mixing, redirecting 
the wave mixing in forward or backward time. This also offers new avenues for flexible control in 
nonlinear optics and potential applications in all-optical computing. 

In an undepleted-pump SHG scheme, we consider the quasimonochromatic waves with car- 
rier frequencies of fundamental wave at u\, and second harmonic wave at u 2 = 2oj\. After repre- 
senting the electric field as £%{z, t) = Ei(z) exp(-iuiit) + c.c = Ai(z) exp[— i(uit — hz)} + c.c., 
where i — 1,2 refers to FW and SH, respectively, and Ai(z) is the slowly varying amplitude (slow- 
varying envelope approximation), the wave equations governing the spatial field envelope E^z) in 
a nonmagnetic nonlinear medium can be current charges expressed as^: 

dA 2 _ UV%deff A 2JAkz 



dz kod 
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-Afe 1 ^ 2 , (1) 



where the phase-mismatching vector Ak = 2k\ — k 2 . And d e ff is the effective nonlinearity, c is 
the speed of light. To consider its time-reversed counterpart, we take the complex conjugate of Eq. 
0] to get 

dA 2 _ j7T iuJ 2deff A *2JAk'z /™ 

~dz~ ~ 6 ~l^ Al 6 • (2) 
Note that Eq. [2] is the same as Eq. [T] except for an extra phase difference (A0 = 2 — 2(j)\ = n) 
and sign-flipped phase-mismatching vector (Ak' = —Ak). The complex conjugate fields E*(z), 
which represents the backwards fields, satisfy Eq. |2j With such exact time-reversal configurations, 
SHG beams can undergo a backward process with respect to its counterpart in forward-time direc- 



tion under a time reversible environment that excluding magnetism and loss. Here the model is 
based on monochromatic wave or waves with simple envelopes, for more complex ones one has to 
time-reverse theirs waveforms in space and time as well, e.g. through time-reversed mirror ^ or 
nonlinear wave mixing^. 

Generally, in the SHG process, SH signals grow along the propagation in a nonlinear optical 
crystal; however, they can also decay due to the familiar phase-matching problem. Here the phase 
plays an important role: it directs the time axis where energy flows. For example, for scheme SHG 
in quasi-phase-matching (QPM) gratings, one can purposely flip the 2nd-order susceptibility (x^) 
of ferroelectric domains adding a phase jump to SH waves in order to prevent SH falling back 
to FW, encouraging SHG conversion. Fundamentally, this originates from an interference effect, 
similar to the linear CPA, however, such interference is nonlinear one due to nonlinear polarization 
at molecular level between SH and FW beams^, hence, the phase difference is crucial here to 
consider the time -reversed process of a SHG as well. 

Figure [TJ shows the scheme for one time-reversed SHG. Here a pair of two identical thin 
BBO (beta barium borate) crystals cut for Type I SHG@ 1064 nm was used (see Supplement). The 
crystal length L is short enough to ensure pump non-depletion (small signal approximation), i. e. 
A±(z) = Ai(0). The corresponding phase-mismatching vectors Ak (Ak') in BBOl (BB02) can 
be independently tuned by the rotating the crystals. The intensity of the SH after BBOl can be 
obtained by direct integrating Eq. [I] to give ^ I 2 = 2dljj-u 2 If/(n 1 n2eoc 2 )L 2 sinc 2 (AkL/2) = 
I norm sinc 2 (AkL/2), where / is the wave intensity, rii is refractive index and e is the vacuum 
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Figure 1: The schematic of SH generation and annihilation. The wave mixing behavior of SHG process 
in the first nonlinear medium is time-reversed in the second (outlined in dashed box). The condition for time- 
reversed symmetry is obtained by adding a it phase shift to FW/SH phase difference and flipping the sign of 
phase-mismatching vector. For a spontaneously grown SHG (starting from zero) in BBOl, its time-reversed 
process leads to "perfect" SH annihilation in BB02 under time reversal symmetry. 
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permittivity. SH is generated when AkL/2 is in the range of (— 7r,7r), and the phase difference 
between the generated SH and FW is A<p = AkL/2. Thus BB01 acts dual-roles as an SH source 
generator and a phase controller. Both the pump FW and generated SH are then directly incident 
into the second nonlinear medium BB02, where time-reversed SHG is constructed. For a sponta- 
neously grown SHG (starting from zero) in BBOl, its time -reversed process points to exactly the 
opposite way: SH signals with the appropriate relative amplitude and phase to the pumping beam 
can be totally annihilated when they incident onto a nonlinear crystal, as will be shown in BB02. 

To gain insight into the nonlinear waves mixing properties along the proposed time-reversed 
structure, we calculate the total output after BB02 by integrating Eq. [T]along the propagating path. 
We refer the normalized SH intensity I norm to unity hereafter, which equals to the output of BBOl 
at phase-matching condition. The total SH output is studied by scanning through phase-matching 
vectors for both crystals, as shown in Fig. [2^. The points of interest are those where SH generated 
in BBOl is totally cancelled out in BB02 leaving only FW out of BB02. In Fig. [2^, where the 
two crystals are simply cascaded, SH vanishes if only AkL = Ak'L = tt/2. However, this is well- 
known that SHG at phase-mismatching condition induces energy flow oscillation between FW and 
SH at the period of twice of the coherent length L c = ix / Ak. The four corners are not the desired 
ones, which are just due to non-SH input and also no SH generated by BB02. Figure[2j) shows the 
typical SH intensity along the propagating path at the conditions of phase-matched, quasi-phase- 
matched, and phase-mismatched conditions (shown in red, green and blue lines, respectively). The 
situation of green line resembles the idea of QPM, which flips the sign of phase-mismatching 
vector when FW and SH are out of phase at L c to gain a continuous SH generation. Note that, 
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Figure 2: The total SH output characteristics and the SH intensity along the propagation path, a, The 

normalized total SH output without a phase shift between the two crystals. The two white dots show two 
situations when SH generated in BBOl is cancelled out in BB02 due to dephasing. b, SH intensity along 
the propagation length in the two crystals in phase-matching and mismatching situations, c, The normalized 
total SH output is greatly modified after introducing an additional it phase shift to A0 in between the two 
crystals. The diagonal dashed line shows total SH outputs equal to none, indicating any SH generated in 
the first nonlinear medium can always be cancelled out in the second, d, SH intensity along the propagation 
length in the two crystals. The SH intensity evolution manifests itself in spatially symmetric pattern, which 
is also an indication of time-reversal of the wave mixing process. White dash dots and line in a and c indicate 
zero intensity. 
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SH's intensity only flows back under the phase mismatching condition, and occurs at twice of the 
coherent length. Along its propagation, SH exhibits symmetry between the two crystals. This 
stimulates the quest of exploring time- symmetry for re-converting SH back to FW at any moment 
despite their phase relations. 

According to previous analysis, time reversal symmetry can be implemented by introducing 
an additional n phase difference before BB02 and reversing the phase mismatching vector in 
BB02. We re-plot the SH output after BB02 in Fig. [2];. As it shows, the SH generated in 
BBOl can always be cancelled out in BB02 whenever Ak' = —Ak is fulfilled. This reveals the 
circumstance where SHG generated by BBOl undergoes annihilation process along a time-reversal 
path in BB02. Figure[2]i illustrates this idea better by showing SH intensity along the propagating 
path across both crystals. With these conditions, the intensity profiles of SH preserves the spatial 
symmetric along the propagation axis through two crystals. It is more convincing to observe such 
symmetry even under the phase-matched condition in Fig. [2ji. In all these circumstances, the phase 
plays a core role in the determining of the direction of the wave mixing process. 

Experimentally, we exam this time-reversal SHG by considering varying the phase mis- 
matching vector to satisfy the condition for total annihilation. We set BBOl for several differ- 
ent conditions where AkL/2 = — 0.3 757T, — 0.2257T, 0, 0.357T, 0.67T. The SH intensity generated 
in BBOl is If = sinc 2 (AkL/2) with the FW/SH phase difference A0 = AkL/2. Then a vr 
phase shift to A<p is added between the two BBO crystals (See Supplement). Figure [3] shows the 
measured SH output from BB02 by scanning the phase mismatching vector. Significant annihila- 
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Ak'L/2 (in unit of tt) Ak'L/2 (in unit of tt) Ak'L/2 (in unit of n) 

Figure 3: Experimental measurements of reversed SHG. The measured total SH intensity with different 
FW/SH inputs into BB02 by varying its phase-mismatching vector Ak'. The dips are the signature of 
CPA, where SH generated in BBOl is significantly absorbed by time-reversed SHG in BB02. (Insets) 
Photographs of the SH spot at the specified conditions. 

tion of SH in BB02 occurs when the time-reversal requirements are met, that is, the dips of the 
SH output curves are located at Ak' = —Ak. The results clearly indicate that for an SHG even 
with low conversion efficiency there always exists a symmetrically reversed SHG process that is 
capable of "perfect absorbing" its SH wave when coherently illuminated by both FW and SH. The 
corresponding attenuation to each dip of the curves is measured to be 11.5 dB, 12.4 dB, 14.5 dB, 
1 1.4 dB and 10.6 dB, yet the SH generated in BBOl is predicted to be totally absorbed in BB02 
theoreticallyE^GII. Such modulation depth is limited by the experimental conditions, e.g. pumping 
intensity, wavefront, coherence. However, strictly complete annihilation cannot be reached should 
quantum mechanism and quantum noise are considered 1 ?^. 

This kind of time-reversal phenomena are rather general and should be expected to occur for 
nonlinear optics with SHG in the up-conversion region, but also another typical down-conversion 
mixing - OPA, which provides an important example in nonlinear optics. Moreover, similar to two- 
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channel CPA, the three waves mixing in an OPA process also adds more complexity by considering 
phases and gain. Compared with reversed SHG, the pumping effect can plays a critical role in 
reversed OPA geometry. Here we extend our study to a reversed two-channel OPA, which is 
qualitatively requires two coherent input beams and the pumping beam. Annihilation of waves can 
only be achieved when the relative phase, amplitude and pumping level conditions are reached, as 
schematically shown in Fig. [4^. Thus, it is not only sensitive to frequency but also to the amplitude 
and phase of the input light. 

We consider non-depleted Type II degenerate OPA of three quasimonochromatic waves with 
carrier frequencies of fundamental wave at ui p , signal and idler at tu s = Ui = co p /2 in BBO crystal. 
By slow-varying envelope approximation, neglecting group velocity mismatch (GVM) and group 
velocity dispersion (GVD), the envelops satisfy the coupled wave equations^ 

dA s _ 2iu 2 s d eff A * iAkz 
~dz~~ ~J^ ApAe ' 

dA i = 2iuJ f d eff A A * e iAkz 

dz /c^c^ 

The subscripts p, s, i stand for pump, signal and idler, respectively. The phase mismatching vector 
Ak = k p — k s — ki is required to be zero for efficient conversion. Hence, the relative phases 
matters the most for now. d e f j is the effective nonlinearity and c is the speed of light. One can 
find the phase-matched parametric gain of signal/idler with respect to pump intensity and total 
phase difference (A0 = 20 p — (f> s — <pi) to be l 16 l 23 lHZ|: $ = \cosh(gz) — iexp(—iA(j))sinh(gz)\ 2 , 
in which g = 2^J(u s uJi) / {n s n^)d e ^A v j c is proportional to the pump amplitude. The results with 
different pump intensities are shown in Fig. |4j). The gain or attenuation is governed by both the 
pump intensity and the total relative phase A0. Hence, the relative phases matters the most when 
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we consider the time reversed OPA. Since A&; = 0, OPA and its time-reversal process is accessible 
in the same condition, an OPA device can behave like an OPA-CPA device^. Theoretically, the 
in-phase signal can be amplified by the factor $, whereas the out-of-phase signal is attenuated by 
the same factor. The dip in the gain curve at A0 = rr/2 is the clear indicator of annihilation of 
Signal/Idler beams when pumping above the threshold of OPA. 




Figure 4: The schematic of OPA and its reversal, a, The schematic of Signal/Idler pair amplification 
and annihilation, b, Behavior of the parametric gain <I> vs. total phase difference Acj) for various pump 
intensities, c, The experimentally measured overall <I> vs. A<p. Gain variety along pump pulses is taken into 
account in data fitting. 

In our reversed-OPA experimental setup, the laser source is a Nd: YAG nanosecond laser sys- 
tem which nominally delivers 4-ns (FWHM) FW pulses at 1064 nm and 3-ns SH (see Supplement). 
The polarization of attenuated FW is rotated to ensure balanced orthogonally polarized Signal/Idler 
input onto the nonlinear medium (BBO). The device is pumped by the SH and its phase is finely 
tuned to adjust the total phase difference. In Fig. |4j:, we show the experimentally measured para- 
metric gain when the phase difference between the three waves is varied. Far from A0 — n/2, the 
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system behaves as an OPA device($ > 1). Note that, at A0 = tt/2, a dip with parametric gain 
$ < 1 can be observed in the gain curves, indicating the Signal/Idler annihilation, when pumping 
above the threshold of OPA. Such a dip is the clear signature of "colored" OPA-CPA 16 , which is 
also the time-reversal process of its counterpart OPA above threshold. In an ideal case, the product 
of the gain factor for the in-phase and out-phase signal, i.e. maximums and minimums on the 
curve, would be a unity. This is a straightforward outcome of time-symmetry rule. However, there 
are plenty of practical issues that would affect this relation. One critical parameter is the pulse 
duration. The actual pulse of pump wave with shorter time duration does not completely overlap 
with the Signal/Idler pulses, meanwhile, the actual gain along the pump pulse varies. Nevertheless, 
these dips still clearly indicate that such reversed-OPA device is capable of attenuating Signal/Idler 
given the coherent illumination with correct phases, which behaves as a "colored" OPA-CPA sys- 
tem. The result demonstrates two operating states of a wave mixing process in its forward and 
time-reversed direction. However, in our experiment, gain depletion regime is never reached. In 
this regime, OPA-CPA is still possible when energy can be drained completely from the pumping 
beamHHl. This requires further studies in the future. 

In conclusion, we have demonstrated time reversed wave mixings for SHG and OPA. This 
enables us to observe the annihilation of coherent beams under time-reversal symmetry by vary- 
ing the relative phase of the incident fields. Time reversed SHG is able to absorb coherent wave 
at second harmonic frequency of its pump. For time reversed OPA, we show that the OPA can 
simultaneously amplify and attenuate coherent signal and idler waves when pumped above thresh- 
old. Our study provides a versatile platform for flexible control in nonlinear optics and potential 
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applications in efficient wavelength conversion, all-optical computing. 
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